ASYMPTOTIC BEHAVIOUR OF THE EINSTEIN- VLASOV 
SYSTEM WITH A POSITIVE COSMOLOGICAL CONSTANT 



HAYOUNG LEE 

Abstract. We study locally spatially homogeneous solutions of the Einstein- 
Vlasov system with a positive cosmological constant. First the global existence 
of solutions of this system and the casual geodesic completeness are shown. 
Then the asymptotic behaviour of solutions in the future time is investigated 
in various aspects. 



1. Introduction 

The Vlasov equation in general describes a collection of collisionless particles. 
Each particle is driven by self-induced fields which are generated by all particles 
together. In gravitational physics, which is the subject we are considering here, the 
fields are described by the Einstein equations. In other words, the Vlasov equation 
is coupled to the Einstein equations by using the energy-momentum tensor. One 
application of the Vlasov equation to this self-gravitating system is cosmology. And 
the particles are in this case galaxies or even clusters of galaxies. 

We study the dynamics of expanding cosmological models in the presence of a 
positive cosmological constant. The simplest cosmological models are those which 
are spatially homogeneous. Spatially homogeneous spacetimes can be classified into 
two types. The conventional terminology is that there are Bianchi models and the 
K ant owski- Sachs models. The models with a three-dimensional group of isome- 
tries G3 acting simply transitively on spacelike hypersurfaces are Bianchi models. 
There are nine types I-IX, depending upon the classification of the structure con- 
stants of the Lie algebra of G3. Those admitting a group of isometries G4 which 
acts on spacelike hypersurfaces but no subgroup G3 which acts transitively on the 
hypersurface are Kantowski-Sachs models. In fact, G3 subgroup acts multiply tran- 
sitively on two-dimensional spherically symmetric surfaces. See |17II19| for general 
information on this classification. 

At this point it is necessary to discuss which solutions of the Einstein- Vlasov 
system should be considered as cosmological solutions. If we take as a cosmological 
spacetime one which admits a compact Cauchy hypersurface, the Bianchi types 
which can occur for a spatially homogeneous cosmological model are only type 
I and IX and also Kantowski-Sachs models. Because of the existence of locally 
spatially homogeneous cosmologies, we take a larger class of spacetimes possessing 
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a compact Cauchy hypersurface so that this allows a much bigger class of Bianchi 
types to be included. Since the Cauchy problem for the Einstein- Vlasov system is 
well-posed, it is enough to define the class of initial data. Here is the definition. 

Definition 1. Let 5^, kij and / be initial data for a Riemannian metric, a second 
fundamental form, and the Vlasov equation, respectively, on a three-dimensional 
manifold M. Then this initial data set (9ij,kij,f) for the Einstein- Vlasov system 
is called locally spatially homogeneous if the naturally associated data set on the 
universal covering M is homogeneous, i.e. invariant under a transitive group action. 

So the spacetimes considered here will be Cauchy developments of locally homo- 
geneous initial data sets on some manifolds. Singer's main theorem in |15j says 

If M is a complete, simply connected Riemannian manifold which is 
locally homogeneous, then M is a homogeneous Riemannian mani- 
fold. 

So one can see that a complete Riemannian manifold is locally homogeneous if 
and only if the universal cover is homogeneous. For Bianchi models the universal 
covering space can be identified with a Lie group G. So the natural choice for G 
in this case is a simply connected three-dimensional Lie group. (For a detailed 
discussion on this we refer to jSJ |Sj ) • 

Without a cosmological constant, it is known that solutions of Bianchi type IX 
and Kantowski-Sachs models cannot expand forever, while the other Bianchi types 
force the volume to be monotone. In fact, if the Bianchi type is IX or the model is 
Kantowski-Sachs, then it has curvature singularities after finite proper time both in 
the past and in the future directions. All other Bianchi types which are expanding 
at some time have a curvature singularity at a finite time in the past and all these 
are future geodesically complete (see JSJE]). 

The detailed picture of the asymptotics of the spatially homogeneous cosmolog- 
ical models without a cosmological constant has not yet been achieved in general. 
However in ^2] for the Einstein- Vlasov system with massless particles the LRS 
(locally rotationally symmetric) reflection-symmetric solutions of Bianchi types I, 
II, III, VIIo and Kantowski-Sachs type have been analysed as far as to give a full 
description of their general behaviour near the singularity and in a phase of unlim- 
ited expansion. Massive particles for LRS reflection-symmetric solutions of Bianchi 
types I, II and III are handled in j!4j . Also in ^H] LRS models of Bianchi type III 
are analysed, which extends and completes the results of ^3] an d In the case 
where the matter model considered is a perfect fluid, much more is known. 

In the presence of a positive cosmological constant, in the existence of inho- 
mogeneous solutions of the Einstein- Vlasov system is proved and future asymptotic 
behaviours are investigated, in the cases of plane and hyperbolic symmetries. The 
spatially homogeneous cases of the results of the plane and hyperbolic symmetries 
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in J^i correspond to LRS models of Bianchi type I and LRS of Bianchi III re- 
spectively. Surveys of results on the Einstein- Vlasov system in general and other 
cosmological models can be found in [11171 151 131 1111 IT2"] . 

When a positive cosmological constant is present, Bianchi IX and Kantowski- 
Sachs models have complicated features. These have both expanding and recol- 
lapsing phases. One can find some results on Kantowski-Sachs models in [3]. In 
3 Kantowski-Sachs model with perfect fluid matter is analysed. In Bianchi 
IX with dust model has been studied. It has been found that there are chaotic 
behaviours between expanding and recollapsing phases. In this paper, in discussing 
expanding cosmological models, we restrict ourselves to Bianchi types except IX. 

At late times Bianchi models except IX with a cosmological constant must ex- 
pand exponentially. This kind of exponential expansion is called inflation. The 
inflationary models of the universe involve a very rapid expansion close to the big 
bang. Also the fact that the expansion of our universe seems to be accelerating 
follows from observations of supernovae of type la. 

The content of the rest of this paper is the following. In Section[3 we discuss the 
global existence of solutions for the Einstein- Vlasov system and the causal geodesic 
completeness of the spacetime towards the future. In Section we investigate 
the asymptotic behaviour of solutions at late times in various aspects. As we will 
see later on, in the case where a positive cosmological constant is present, the 
situations of the locally spatially homogeneous cosmological models have a unified 
picture in their expanding dynamics. We observe the future asymptotic behaviours 
of the mean curvature, the metric, the momenta of particles along the characteristic 
curves as well as the generalized Kasner exponents and the deceleration parameter. 
Also we analyse the energy-momentum tensor in an orthonormal frame on the 
hypersurfaces. 

Here is the formulation of the Einstein- Vlasov system with a positive cosmolog- 
ical constant. Let G be a simply connected three-dimensional Lie group and {e{\, 
a left invariant frame and {e 1 }, the dual coframe. Consider the spacetime as a 
manifold GxJ, where I is an open interval and the spacetime metric of our model 
has the form 



The initial value problem for the Einstein- Vlasov system is investigated in the case 
of this special form of the metric and the distribution function / depends only on t 
and v l , where v l are spatial components of the momentum in the frame e.;. Initial 
data will be given on the hypersurface G x {0}. 
Now the constraints are 



ds 2 = -dt 2 + g i3 {t)e l ® e j . 



(1.1) 



R - + = 16^T 00 + 2A, 



(1.2) 
(1.3) 
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The evolution equations are 

d t9ij = -2k l3 , (1.4) 
dthj = R l3 + (k lm g lm )k tJ - 2k a k) - totTy (1.5) 
- 47rT 00 gy + ^{T lrn g lm )g l:j - kg^. 

These equations are written using frame components. Here fey is the second fun- 
damental form, R is the Ricci scalar curvature and Rij is the Ricci tensor of the 
three-dimensional metric. And A is a positive cosmological constant. 
Here are the components of the energy-momentum tensor ; 

Too(t) = J f(t,v)(l + g rs v r v s )y 2 (detg)^ 2 dv (1.6) 
r«(t) = J f(t,v)vi(detg)^ 2 dv (1.7) 

T ij (t) = J f(t,v)v tVj (l + g rs v r v s )- 1/2 (^tg) 1/2 dv. (1.8) 

Here u := (v 1 , v 2 , v 3 ) and dv := dw 1 c£i> 2 cfo> 3 . 
The Vlasov equation is 

d t f + {2ky - (1 +g rs v r v s )- 1/2 J t mn v m v n }d v ,f = (1.9) 
Here the Ricci rotation coefficients r f mn are defined as 

linn = 2^ (~ C l nk g m i + C [ km g n i + C l mn gki) 

where Cj fe are the structure constants of the Lie algebra of G. 

In general the evolution equations are partial differential equations. However 
since the spacetime is locally spatially homogeneous in this paper, this gives a big 
simplification. I.e., dt is actually d/dt in the evolution equations. Yet the Vlasov 
equation remains as a partial differential equation. We abuse the notation in this 
way to follow the notation in the literature. 

2. The global existence of solutions and 
the future geodesic completeness 

In the absence of a cosmological constant the existence of solutions for the 
Einstein- Vlasov system has been proved in [5] by A. D. Rendall under some condi- 
tions. First under the following assumptions on the initial data : 

(1) gij(O), fcy(O) and f{0,v) are initial data for equations i|1.4l> . I1.5fl and (|1.9|) 
which have Bianchi symmetry and satisfy the constraints (11.21) and (|1.3|) . 

(2) f(0,v) is a C 1 function of compact support, 

the local existence of solutions has been proved. The argument to achieve the global 
existence is that, if \g\, (det g)" 1 , \k\, ||/||oo and the diameter of the support of / 
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are bounded on the finite time interval, [0, T) where the solutions exist locally, then 
T = oo. Here g and k are matrices with components gij and kij respectively. And 
\g\ denotes the maximum modulus of any components g^ for all t, with the similar 
definition for \k\. Another criterion for the global existence, which is a stronger 
version, is the following. If the mean curvature is bounded on [0, T), then \g\, 
(detg) -1 , \k\, \\f\\oo an d the diameter of the support / are also bounded in [0, T). 
In another words, the boundedness on the trace k implies T = oo, i.e., solutions 
of the system exist globally. Here the existence of solutions means the existence of 
solutions and the uniqueness. 

In the case of the presence of the cosmological constant, to obtain the local 
existence of solutions, these results from can be directly applied with a minor 
change. So it is enough to obtain the boundedness of the mean curvature to achieve 
the existence of the global solutions. 

Here is the statement of the global existence of solutions for the Einstein- Vlasov 
system with a cosmological constant. 

Theorem 1. Let (g^ -(0), fcy(0), /(0, v)) be an initial data set for the evolution 
equations \l-4\l , ^1-5}) and the Vlasov equation fl.yp which has Bianchi symmetry 
and satisfies the constraints and \1.3\) . Also let f(0, v) be a nonnegative ex- 

junction with compact support. Then there exists a unique C 1 solution (g^ , kij, /) 
of the Einstein- Vlasov system, for all time t. 

Proof. <|1.4|) - (|1.9|) imply a homogeneous first order ordinary differential system 
for l|1.2(l and (|1.3|l . To be more precise, define 

A :=R — {kijk ij ) + (k tj g lj ) 2 - 16irT QQ - 2 A, 
Ai := V i k ij + 8nT 0j . 

Then after a lengthy calculation we obtain 

d t A = 2{k ij g»)A-2 1 \ j g»A l , 
d t Ai = (k lm g lm )A l + 2k l l A l . 

This means that if the initial data satisfy the constraints l|1.2|l and l|1.3|l then so 
does the solution of the evolution equations (|1.4(l and (|1.5(l with energy-momentum 
tensors (f 1 . f>|) - (f 1 . S|) and the Vlasov equation (|1.9(l . As we have mentioned before, by 
the same argument in [S] (see pp. 86-89), the solution is uniquely determined by 
the initial data on [0, T), for some time T. So we need to check the boundedness 
of the mean curvature to achieve the global solution. This follows from . 2|> in 
Proposition^ which we state and prove in the next section where it is mainly used. 
□ 
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The next result asserts the geodesic completeness of locally spatially homoge- 
neous spacetimes for the Einstein- Vlasov system with a positive cosmological con- 
stant. 

Theorem 2. Suppose the hypotheses of Theorem^ hold. Then the spacetime is 
future complete. 

We postpone the proof of this theorem to the next section. 

3. The future asymptotic behaviour of solutions 

Bianchi models except type IX with a cosmological constant correspond to ex- 
ponential inflation. Also the spacetime becomes isotropic at late times. In other 
words, all locally spatially homogeneous solutions of the Einstein- Vlasov system 
with a positive cosmological constant are approximated by the de Sitter solution as 
time goes to infinity. The idea that the de Sitter solution should act the attractor 
for solutions of the Einstein equations with a cosmological constant is called the 
cosmic no hair conjecture. Wald's result fjlSp can be considered as a proof of this 
conjecture in the homogeneous case. In |18| the form of the energy-momentum 
tensor is not specified in detail. It satisfies only the dominant and strong energy 
conditions given respectively by 

(1) T a! 3V a w^ > where v a and w@ are any two future pointing timelike vectors, 

(2) R a f}V a v^ > for any timelike vector v a . 

Here Greek indices take the values 0, 1, 2, 3. 

In this section, we investigate the asymptotic behaviour of solutions of the 
Einstein- Vlasov system in the future time. The basic ideas of the following propo- 
sition are from [T^J. Here we carry out the computation carefully so that the error 
terms are explicitly determined for the future reference. In general throughout the 
paper, a constant C is positive and may change line by line. 

Proposition 1. Let cry be the trace free part of the second fundamental form kij 
such that 

hj = ^{kimg lm )g l j + o-ij. (3.1) 

Then we have 

(hjg ij ) = -(3A) 1 / 2 + 0{e- 2 ^) (3.2) 
(a^) = 0(e~^) (3.3) 

where 7 2 = A/3. 

Proof. From evolution equations (|1.4|) and (|1.5|l . we have 

dt(k ijg ij ) =R+ (%9 U ) 2 + MT t] g tJ ) ~ 12^T 00 - 3A. (3.4) 
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Using the constraint 1)1. 2[l . ()3.4)l becomes 

9 t (% 5 y ) = (ktj k ij ) + \-n,A,' J ) + 4ttT 0Q - A. (3.5) 

By the dominant and strong energy conditions and using 1)3. we obtain from 
(|3.5|) the following : 

d t {k ijg ij ) > (kij k^ ' ) - A > i (fcy s « ) 2 - A. (3.6) 

Also we rewrite the constraint (11.2)1 using (|3.1|) as 

\(hi9 ij ? = -\R + + SttTqo + A. (3.7) 

In j Wald has proved that in all Bianchi models except type IX, the Ricci scalar 
curvature is zero or negative. So from l|3.7|l . we can obtain that the last quantity 
in (|3.(i|l is nonnegative, i.e. 

dt(k tl g^) > ^(k ij9 ij ) 2 - A > 0. (3.8) 

From the last inequality, one can see that if the universe is initially expanding, i.e., 
the mean curvature is nonzero initially which will be the case here, then the mean 
curvature never becomes zero and this means that the universe expands forever. 
Now the first case for the solutions from (|3.8(l is by taking the equal sign in the last 
inequality : 

{k ij9 ») = -(3A) 1 / 2 (3.9) 

for all time t. We choose the negative sign as the convention throughout the paper. 
The other case is by solving the following ordinary differential inequality : 

dt(k ijg v) < 1 
(%^) 2 -3A " 3' 

I.e., integrating both side of the following 

/ (fc^)-(3A)V 2 \ 
nN^) + (3A)V 2 J^ 

we obtain 

> (%g y ) + (3A) 1 / 2 > -Ce~ 2 ^ 
where C is a positive constant. Therefore combining with (|3.9|) . we get 

(^■^) = -(3A) 1 /2 + o (e -27t ) . (3 . 10 ) 

Here while we use the notation O(-), we lose a piece of information from l|3.8(l that 
the mean curvature is non-decreasing. So we want to point out that in 1)3.10)1 the 
error 0(e~ 2jt ) is possibly zero or negative error. 

Now to prove ()3.3)) . we use 1)3.7)1 once more. Then we get 

i(^^")<i(fc^") 2 -A. 
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Using [|3.2jl . one can see that 

a 

To identify cr^- from (|3.3|) which is necessary to investigate the asymptotic be- 
haviour of gij , let us use some properties of linear algebra which can be found in 
0. Let Ai and A 2 be n x n symmetric matrices with A\ positive definite. Also let 
A be a n x n matrix. Define a norm of a matrix by 

||A|| := sup{||Ar||/||a;|| :i/0,i£ M™}. 

Also define a relative norm by 

||A 2 || Al := sup{p 2 x||/||A ia; || :i/0,i£ M™}. 

Then from these definitions, one can sec that 

114.11 < H^alUJIAxH (3.11) 

and also 

Malk < {tr(A^A 2 A^A 2 )) 1/2 . (3.12) 
Here is a lemma which gives a relation between <xy and g^ . 

Lemma 1. Let \\g(t)\\, \\k(t)\\ and \\o~(t)\\ denote the norms of the matrices with 
entries g^, kij and a^, respectively. Then 

\\o(t)\\<Ce-^\\g{t)\\ 

where C is a constant. 

PROOF. By (|3.12|l and (|3.3|) we have 

So the lemma follows by l|3.11|l . □ 

The next lemma gives some information on and g % K It does not say the 
full information about the asymptotic behaviours of gij and g lJ . It is a step in an 
argument which leads to stronger estimates on these quantities. 



Lemma 2. 



|e- 2 ^|<C, (3.13) 
|e 2 TV J '|<C, (3.14) 



for all t > 0. Here C is a constant. 
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PROOF. Let g^ := e 2yt gij- Then using (|3.1jl and (|3.2|) . we get 

2 

= 0{e-^ t )g ij -2e-^ t a ij . (3.15) 

Now let us use the norms again. Let be a norm of the matrix with entries 

g~ij. Then using (|3.11|) we get 

||ff(*)ll< ||5(0)||+ f Ce- 2 -* s (\\g( S )\\ + \\a( S )\\)d S 
Jo 

< ||.g(0)|| + f Ce-^(\\g(s)\\ + [k(*)|| sW ||fl(*)||) da. 
Jo 

Note that 

So with (|3.3|) we have 

H5(*)ll< 115(0)11+ f Ce-^\\g{s)\\ds. 



By Gronwall's inequality, we obtain 

\\g(t)\\ < ||5(0)||exp 



Ce~"> s ds 



< C 



where C is a constant. Therefore gtj is bounded by a constant for all t > 0. (|3.14|) 
can be proved by the same argument. So we omit it. □ 
Combining Lemmas ^ an(1 [21 wc obtain the following proposition. 

Proposition 2. 

/or Zarge t. 

The Bianchi models except type IX with a cosmological constant are inflationary 
cosmological models. I.e., the metric g^ increases exponentially which we have 
already seen in Lemma [21 The next proposition gives us detailed information on 
the error term of g^ which allows us to construct g lJ . I.e. in order to find g lJ , 
determining the form of the error term in g^ is crucial. The first main point of 
this proposition is to find g^ with the error term, using Propositions ^ and [21 And 
then we want to obtain g lJ explicitly which will play an important role in the rest 
of the paper. 

Proposition 3. 

gij {t)=e 2 ^(g ij + 0(e-^)) (3.16) 
g*3(jb) = e- 2 ^^ + 0{e-^)) (3.17) 

for large t. Here Gij and Q lJ are independent of t. 
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Proof. Consider (|3.15(l again with Lemma[21and Proposition Then we have 
dtgij = 0(e- 2 ^)g l3 - 2e- 2 ^a l3 
= 0(e~~ ft ). 

Since dtgij is decaying exponentially, there exists a limit, say Qij , of g~ij as t goes 
to infinity. Then this gives 

9ij = Qij + 0(e-^) 

i.e. 

g ij = e ^t{g ij + {e-^)). 

Here the lower order term of gij is of an exponential form, i.e. C(e 7 *) so that it 
is combined with the leading order term, which makes it possible to compute g v 
explicitly. So g^ is 

g v = e - 2 ~< t (g ij + 0(e-~< t )y 1 

for large time t. Note that the last equality is due to Taylor expansion. □ 
Now let us consider the generalized Kasner exponents. Let A^ be the eigenvalues 
of with respect to gy, i.e., the solutions of det(fc*- — AS*) — 0. Define the 
generalized Kasner exponents by 

_ Ai Xi 
Pl ~ (W m )' 

The name comes from the fact that in the special case of the Kasner solutions these 
are the Kasner exponents. Note that while the Kasner exponents are constants, the 
generalized Kasner exponents are in general functions of t. The generalized Kasner 
exponents always satisfy the first of the two Kasner relations, but in general do not 
satisfy the second, where these two Kasner relations are 

I> = 1, (3.18) 

i 

E^) 2 = L ( 3 - 19 ) 

i 

The following proposition says that all generalized Kasner exponents Pi go to ^ as 
t goes to infinity. This means that the spacctime isotropizes at late times. 

Proposition 4. 

p. l {t)= l - + 0{e-~< t ) 

for large t. 

Proof. First note that by l|3.1|l Xi are also the solutions of 
det {a i j -[X-^(k lm g lm )}5^)=0. 
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So the eigenvalues of aij with respect to gij are 

A, :=\i-\{k lm g lm ). 

Also note that Ei(A;) 2 = oi m a Xm ■ Then by we get 

Aj = C(e-T*). 

Therefore with ll-i.2t we obtain 

Pi~\ = i„ K lm =0{e-^). 
3 s{ki m g lm ) 

□ 

Before we look at the momenta of the distribution function, there is another 
quantity to be mentioned regarding expanding models, which is the deceleration 
parameter, say q. This deceleration parameter is related to the mean curvature, as 
followings 

d t (k ij g i i) = -(l + q)(k ij g ij ) 2 . 

In the inflationary models, the deceleration parameter is negative. Especially ex- 
ponential expansions lead to the value —1 for the deceleration parameter. Here is 
the statement. 

Proposition 5. 

q = -l + 0(e- 27t ) 

for large t. 

Proof. It follows by Proposition^ D 
Now we will investigate the behaviour of the momenta of the distribution function 

/ along the characteristics where / is a constant. From the Vlasov equation (|1.9|) 

we define the characteristic curve V l {t) by 

— = - (1 + g r sV r V s )- 1/2 l l mn V m V n (3.20) 

for each V l (t) = v % given i. The characteristics Vl, rather than V 1 , have a bit 
simpler form, which makes investigating the behaviour of the momenta easier. So 
here Vi{t) satisfies 

^ = -(1 + 9rsV r V s r 1/2 lL n V P V q gP m g^g l3 (3.21) 

for each Vi(i) = Vi given i. For the rest of the paper, the capital V % and Vi indicate 
that v % and Vi are parameterized by the coordinate time t, respectively. 

Theorem 3. Consider all Bianchi types except type IX. Then Vi(t) from \S. 21\) 
converges to a constant along the characteristics, as t goes to infinity. I.e., V t (t) — 
e~ 27 *(C + C(e~ 7 *)) where C is a constant, for large time t. 
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This theorem says that the particles in the cosmological models with a positive 
cosmological constant are expanding with a momentum whose components v l decay 
exponentially for future times. 

Before proving Theorem |3 we require some lemmas. 
Lemma 3. 

where V is a constant and for large time t. 

Proof. First note that by Propositions and we have 

a ij = C(e~ 37 *). 

Let <r y := e 37 *cr u . Then is bounded by a constant for all t > 0. Since Q tJ in 
Proposition |21 is positive definite and time independent, there exists a constant C, 
independent of time, such that 

a ■ (V'\,\ r 

Then combining this and l|3.2[l . we obtain 
d_ 



= l(ki m g lm )g lJ v t v 3 + 2^V t V 

< ( - 2 7 + Oie-^g^VtVj + Ce'^G^ViVj 

< ( - 2 7 + Ce-^g^V.Vj. (3.22) 

Here to get the first equal sign, H3.21JI is used. Yet the terms involved with (|3.21(l 
vanish due to the antisymmetric property of ~f mn combining with g 1 ^ . Now consider 
V := e^g^ViVj. Then one can see from (l3~22j) that 

^ = 0{e~^)V. (3.23) 
at 

So there exists a limit of V, say V as t goes to infinity. This gives that 

gi'ViVj = e- 2 ^(V + 0(e'^)). (3.24) 



□ 



Lemma 4. 

\Vi\ < c + o( e - 7t ) 

where C is a constant, and for all i and for large time t. 

Proof. Since Q % i is positive definite, there exists a constant C such that 

\Vi\ 2 < cg^ViVj. 

With (|3.17(l . Lemma |21 implies that 



Asymptotic Behaviour of the Einstein- Vlasov system 



13 



So the lemma follows. □ 
Proof of Theorem[31 Note that in Bianchi type I, since all structure constants 

are zero, also the Ricci rotation coefficients 7„ m are zero. So from (|3.21|> it is clear 

that Vi(t) — vi for all t, where V{ = Vi(t) for a given time t. 

More generally Lemma 0] tells us that all Vi are bounded by a constant when t 

goes to infinity. However this allows oscillating behaviours. So to rule out these 

cases, we need to go back to look at Using l|3.16[l and (|3.17() . we have 

dV 



dt 



<Ce- 2 *\V p \\V q \ 



Here the right hand side is a summation for some p and q. So by Lemma 0] we 
obtain 

dt 

where C is a constant. This leads to the conclusion that when t goes to the infinity, 
goes to zero exponentially, i.e. Vi goes to a constant. Also by l|3.17[) . it is clear 
that V l {t) ^e~ 2 ^(C + 0{e~^)). □ 
Now let us prove the geodesic completeness we have postponed in the previous 
section. 

Proof of Theorem [3 The geodesic equations for a metric of the form 
imply that along the geodesies the variables t, v l and v° satisfy the following system 
of differential equations : 

t =V °' (3 " 25) 

dv ° U i 3 

— = 2kyv° - 7 ^"V 

where r is an affine parameter. For a particle with rest mass m moving forward in 
time, v° can be expressed by the remaining variables, 

v° = (m 2 + PijuV) 1 / 2 . (3.26) 

The geodesic completeness is decided by looking at the relation between t and 
the affine parameter r, along any future directed causal geodesic. This relation is 
clear from (|3.25|) and (|3.26|) . I.e., it is given by 

^ = (rn 2 +. 9y W) 1 /2. 

To control this, we need to control gijV % v^ as a function of the coordinate time t. 
Consider first the case of a timelike geodesic. I.e., m > 0. Then satisfy 

dV i 

^ = 2k)v' - (m 2 + j#vr 1/2 7Lrr. 
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So we can use the proof of Lemma El and we get gijV l V : > is bounded. Hence this 
gives that 

^ = (m 2 +g^V i V j )- 1 / 2 >C (3.27) 

where C is a positive constant. Therefore when r is recovered by integrating i|3.27[) . 
the integral of the right hand side diverges as t goes to infinity. 

Now consider a null geodesic, i.e., m = 0. Then in this case, V l (t) satisfy 

dV i 

— = 2k$V* - (g l] V^)- 1 ' 2 1 l nn V m V n . 
So the proof of Lemma|3]is still applied. Again gijV % V^ is bounded and 



J = faVV)- 1 / 9 > c 



where C is a positive constant. Therefore as t goes to infinity so does r. □ 
The rest of the paper, we analyse the energy-momentum tensor. 

Proposition 6. Let us consider the energy-momentum tensor in an orthonormal 
frame, {e{\. The energy-momentum tensor is described by 



p(t)= J f(t,v)(l + \v\ 2 ) 1/2 dv (3.28) 
Ji(t) = J f(t,v)vidv (3.29) 
Sij(t) = I f(t,v)viVj(l + Ivfy^dv (3.30) 



where p := Too *s the energy density, Ji := To,; the components of the current density 
and Sij := Tij are the spatial components of the energy-momentum tensor. Here 
the hats indicate that objects are written in the orthonormal frame. Furthermore 
v := #2, $3) and dv — diidv2dv3 . 

Then pit), Jj(t) and Sij(t) tend to zero as t goes to infinity. More precisely, 



p(t) 


= 0{e 


-3 7t)) 


Ji(t) 


= 0{e 


-4 7t)) 




= 0{e 


-5 7t)) 



for large time t. Furthermore 

Ji(t) 
Pit) 



C(e~ 7t ), (3.31) 



= 0(e- 2 ^), (3.32) 



for large t. 
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Proof. Combining with the fact that f(0,v) has compact support on v, The- 
orem OH implies that Vi(t) is uniformly bounded in large time t. I.e., there exists a 
constant C such that 

/(*,») =0, if H>C* (3.33) 
for all large t. By H3.16J) we have 

Vi = 0(e-T*)«j. (3.34) 

So using (|3~33|) and we get 

= / /(t,fi)(l + |fi| a ) 1/a dfi. 

J\vi\<Ce— r* 

Note that since /(t, z)) is a constant along the characteristics, 

\f(t,v)\ < H/oll :=sup{|/(0,T))| :V0}. (3.35) 

So we obtain 



p(t)<C f(t,v)dv 

J\v i \<Ce-"' t 

< C||/ ||e- 3 ^. 



Also by l|3~33|l - J3~5S|) we have 



Ji(t) = / f{t,v)vidv 

J|Ci|<Ce-Tt 



Similarly we get 



< Ce" 7 * / f(t,v)dv 

J\v i \<Ce—' t 

<C||/o||e- 47 '. 



Sy(t)= / /(;,i)^(i + H 2 r 1/2 d{) 

J|«i|<Ce-^* 



< Ce^ 27t / f(t,v)dv 

J\ii i \<Ce—' t 

<C||/o|| e ~ 57t . 

Now let us estimate the ratios Jj/p and Sij/p. Again by (|3.33(l and (|3.34(1 . first we 
get 

Ji(t) J f(t,v)vidv 



p(t) ff(t,v)(l + \v\ 2 ) 1/2 di 
< Ce 



J f(t,v)(l + \v\ 2 ) 1 / 2 dv 
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And similarly we obtain 

Sjjjt) = J f{t,v)v l v ] (\ + \v\ 2 )- 1 ' 2 dv 
p(t) ff(t,v)(l + \v\*y/*dv 

□ 

In this proposition since all components of the energy momentum tensor in an 
orthonormal frame go to zero exponentially as t goes to infinity, one can say that 
in a certain sense solutions of Einstein- Vlasov system with a positive cosmological 
constant are approximated by vacuum Einstein solutions. In a more detailed level 
(I3.31|l and l|3.32|l resemble the non-tilted dust-like solutions (Ji(t) and Sij(t) are 
identically zero). 
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